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In information theory, the link between continuous information and discrete information is es- 
tablished through well-known sampling theorems. Sampling theory explains, for example, how 
frequency-filtered music signals are reconstructible perfectly from discrete samples. In this Letter, 
sampling theory is generalized to pseudo-Riemannian manifolds. This provides a new set of mathe- 
matical tools for the study of space-time at the Planck scale: theories formulated on a differentiable 
space-time manifold can be completely equivalent to lattice theories. There is a close connection to 
generalized uncertainty relations which have appeared in string theory and other studies of quantum 
gravity. 
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It is generally assumed that the notion of distance loses 
operational meaning at the Planck scale, lp ss lCP 35 m 
(assuming 3+1 dimensions), due to the combined effects 
of general relativity and quantum theory. Namely, if one 
tried to resolve a spatial region with an uncertainty of 
less than a Planck length, then the corresponding mo- 
mentum uncertainty should randomly curve and thereby 
significantly disturb the very region in space that was 
meant to be resolved. It is expected, therefore, that the 
existence of a smallest possible length, area or volume, 
at the Planck scale or above, plays a central role in the 
yet-to-be-found theory of quantum gravity. 

In the literature, no consensus has been reached as 
to whether this implies that space-time is discrete. On 
the one hand, quantization literally means discretization, 
and space-time discreteness is indeed naturally accommo- 
dated within the functional analytic framework of quan- 
tum theory, see, e.g., 0. Also, most interacting quan- 
tum field theories are mathematically well-defined only 
on lattices. On the other hand, within the mathematical 
framework of general relativity, space-time is naturally 
described as a differentiable manifold and deep princi- 
ples such as local Lorentz invariance would appear to be 
violated if space-time were discrete. 

There is the possibility that the cardinality of space- 
time is between discrete and continuous, but it is strongly 
restricted by results of Godel and Cohen. In 0, they 
proved that both are consistent with conventional (ZF) 
set theory: to adopt an axiom claiming the existence of 
sets with intermediate cardinality or to adopt an axiom 
claiming their non-existence. Therefore, it is not possi- 
ble to explicitly construct any set of cardinality between 
discrete and continuous infinity from the axioms of con- 
ventional set theory. If space-time is describable as a 
set and if this set is of intermediate cardinality, then its 
description cannot be constructive and requires mathe- 
matics beyond conventional set theory. 

In this Letter, we consider a simpler possibility. In a 
concrete sense, space-time could be simultaneously dis- 
crete and continuous. Namely, in the simplest case, phys- 



ical fields could be differentiable functions which pos- 
sess merely a finite density of degrees of freedom. If 
such a field's amplitude is sampled at discrete points 
of the space-time manifold then the field's amplitude at 
all points in the manifold are reconstructible from those 
samples - if the sample points are spaced densely enough. 
The minimum average sample density which allows the 
reconstruction of fields could be, for example, on the or- 
der of the Planck scale. All physical entities such as 
fields, Lagrangians and actions can then be written either 
as living on a differentiable manifold, thereby display- 
ing external symmetries or, equivalently, as living on any 
one of the sampling lattices of sufficiently small average 
spacing, thereby displaying ultraviolet finiteness. Such 
theories need not break any symmetries of the manifold 
because among all sufficiently tightly spaced lattices no 
particular lattice is preferred. 

In the information theory community, the mathematics 
of classes of functions which can be reconstructed from 
discrete samples is well-known, namely as sampling the- 
ory. Shannon, in his seminal work [3j, introduced sam- 
pling theory as the link between continuous information 
and discrete information. Our aim here is to extend 
this link between discrete and continuous information to 
curved space-times. 

The use of sampling theory in quantum mechanics was 
first suggested in jj] , where a close connection was shown 
to generalized uncertainty relations that had appeared in 
studies of quantum gravity and string theory, see . In 
the simplest case, these uncertainty relations are of the 
form AxAp > -|(1 + /3(Ap) 2 ) and imply the existence of 
a finite minimum uncertainty in position Ai m „ = hy/]3, 
as is easily verified. In Q , it was shown that a finite lower 
bound to the position uncertainty implies that the wave 
functions possess the sampling property, i.e. that they 
can be reconstructed everywhere from discrete samples if 
those samples are taken at a spacing that is at least as 
small as the minimum position uncertainty. In , it was 
suggested that the freedom of choice of sampling lattice 
may be related to gauge symmetries. 
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In recent proceedings, see j7|, sampling theory for 
physical theories on curved spaces was outlined. Building 
on , this Letter introduces sampling theory on Rieman- 
nian and pseudo-Ricmannian manifolds. We will obtain a 
covariant information density cutoff together with a new 
sampling theoretic principle that is consistent with the 
Lorentz contraction of sampling lattices. 

The basic sampling theorem goes back to Cauchy 
in the early 19th century, see [8j|. Consider the set 
of square integrable functions / whose frequency con- 
tent is bounded by to m ax, i.e., which can be written as 
f( x ) — J_J* ] {io)e wx dio . These / are called bandlim- 
ited functions with bandwidth uj m ax- If the amplitudes 
{f(x n )} of such a function are known at equidistantly- 
spaced discrete points {x n } whose spacing is ir/uJmax, 
then the function's amplitudes f(x) can be reconstructed 
for all x. The reconstruction formula is: 

f( X ) = jr f(x n ) si l ( *~^" )aw] (i) 



n— — oo 



(x - x n )uj n 



This sampling theorem is in ubiquitous use, e.g., in dig- 
ital audio and video as well as in scientific data tak- 
ing. Sampling theory, see Q, studies generalizations 
of the theorem for various classes of functions, for non- 
equidistant sampling, for multi-variable functions and it 
investigates the stability of the reconstruction in the pres- 
ence of noise. 

Following 0, we now define a general framework for 
sampling on Ricmannian manifolds. The key assumption 
in the basic sampling theorem is a frequency cutoff, more 
precisely a cutoff of the spectrum of the self-adjoint differ- 
ential operator —id/dx. On a multi-dimensional curved 
space, a covariant analog of a bandlimit is the cutoff of 
the spectrum of a scalar self-adjoint differential opera- 
tor. As an explicit example one may choose the Laplace- 
Beltrami operator A = \g\^ 1 ^ 2 dig l ^\g\ 1 / 2 dj where |(/| is 
the determinant of the metric tensor g. 

Consider then the Hilbert space H of square integrable, 
say scalar, functions over the manifold and the dense 
domain T> C TL on which the considered operator, say 
the Laplacian, is essentially self-adjoint. Using sloppy 
but convenient terminology we will speak of all points 
of the spectrum as eigenvalues, A, with corresponding 
"eigenvectors" |A), keeping in mind that the manifold 
will generally be noncompact and its spectrum therefore 
not discrete. We use the notation | ) in analogy to Dirac's 
bra-ket notation, but with round brackets to distinguish 
from quantum states. The \(f>) that we consider here could 
be, for example, the scalar fields that are being integrated 
over in a quantum field theoretical path integral. 

The operator —A is positive and its spectrum is an 
invariant of the manifold. A spatially covariant "band- 
width" cutoff in nature then means that physical fields 
are elements of V p h — P.T>, where P projects onto the 
subspace of T> which is spanned by the eigenspaces of — A 



with eigenvalues smaller than some fixed maximum value 
Amaz, which could be, e.g., on the order of 1/1%. 

For example, in quantum field theoretical actions this 
type of cutoff arises if —A is the lowest order term 
in a power series in —A whose radius of convergence 
is finite, say 1/lp- Examples are the geometric series 

Iff ZZiHp&r^ and l P 2 En=i( 7 lp<f>*^<f>) n . Such 
series correspond to Planck-scale modified dispersion re- 
lations, a concept that has recently attracted consider- 
able attention in the context of the transplanckian prob- 
lems in black hole radiation and inflationary fluctuations, 
see, e.g., Interestingly, also the Dirac-Born-Infeld 

action may be viewed as providing a minimum length 
cutoff through this mechanism, namely when expanding 
the square root in its action as power series with finite 
radius of convergence. See, in particular, 

If, by this or another mechanism, the yet-to-be-found 
theory of quantum gravity does yield a bandwidth cutoff, 
how do the fields in the physical domain T> p h acquire the 
sampling property? For simplicity, assume that one chart 
covers the A^-dimensional manifold. The coordinates Xj, 
for j = 1, N act as multiplication operators that map 
scalar functions to scalar functions: Xj : 4>{x) — > Xj(j)[x). 
On their domain within the Hilbert space Ti. these oper- 
ators are essentially self-adjoint, with an "Hilbert basis" 
of non-normalizable joint eigenvectors {\x)} with contin- 
uum normalization 1 = fd N x Ig] 1 ^ 2 \x)(x\. We have 
(x\<f>) = 4>(x). Since A, being a differential operator, 
cannot commute with the position operators Xj we ob- 
tain the situation P\x) ^ \x). Thus, on the restricted 
domain T> p h, the multiplication operators xj are merely 
symmetric but not self-adjoint (intuitively, for lack of 
eigenvectors) . Correspondingly, the uncertainty relations 
are modified, similar to the toy cases discussed in 

Consider now a physical field, i.e., a vector \<f) £ T> p f L . 
Assume that the field's amplitudes 4>(x n ) = (x n \(j)) are 
known at discrete points {x n } of the manifold. While 
all position eigenvectors |a;) are needed to span 7i, suffi- 
ciently dense discrete subsets {|£„)} of the set of vectors 
{P|a;)} can span D p h- A field's coefficients {(f>(x n )} then 
fully determine the Hilbert space vector \<jj) £ D p h and 
they determine, therefore, also (x\<f>) for all x. Namely, 
defining K n \ = (x n \X), the set of sampling points {x n } 
is sufficiently dense for reconstruction iff K is invertible. 
To see this, insert the resolution of the identity in terms 
of the eigenbasis {|A)} of —A into (x\(j>): 
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(x\\)(\\4>) d\ 



(2) 



We use the combined sum and integral notation since the 
spectrum of — A may be discrete and/or continuous (the 
manifold M. need not be compact) and it is understood 
that eigenvalues can be degenerate. With K invertible, 
one obtains {X\4>) — ^ n K^ l n 4>{x n ) which, when substi- 
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tuted back into Eq[21 yields 

4>(x) = ^G(x,x n ) <f>(x n ), (3) 

n 

with the reconstruction kernel: 

G(x,x n ) = ^ {x\X)K^ n dX (4) 
\\\<\ ma * 

In conventional applications of information theory, the 
sample points {x n } are required to be dense enough to 
allow stable reconstruction in the presence of noise: func- 
tions reconstructed from small samples must have small 
norm, in the sense that there exists a C > such that 

W) < C^l^")! 2 for all |0) G I?^. (5) 

n 

The minimum sample density for stable reconstruction 
in flat Euclidean space equals the bandwidth volume in 
Fourier space, up to a constant, as was shown by H.J. 
Landau in 13]. It is also the density of degrees of free- 
dom, defined as the dimension of the space of bandlimitcd 
functions which possess essential support in a given vol- 
ume. Reconstruction stability is highly nontrivial already 
for classical signals in flat space: a bandwidth cutoff does 
not prevent bandlimited functions from oscillating arbi- 
trarily fast in an arbitrarily large region. These super- 
oscillations recruit degrees of freedom from outside the 
considered region at the expense of reconstruction stabil- 
ity. As was shown in [l4|, superoscillations in quantum 
mechanical wave functions produce effects that raise ther- 
modynamic and measurement theoretic issues. It should 
be most interesting, therefore, to explore for physical 
fields in curved space-time the role of quantum fluctua- 
tions as noise and to use our new approach to generalize 
Landau's theorem to Riemannian manifolds. For fixed 
noise, the density of degrees of freedom then yields the 
maximum Shannon information density as usual, see Q. 

We close the case of Riemannian manifolds with several 
remarks. In the case of one dimension and equidistant 
samples, Eg s I2I4I reduce to the basic sampling theorem of 
Eqffl Another simple case is that of compact manifolds. 
Their Laplacian possesses a discrete spectrum whose cut- 
off renders T> p h finite-dimensional. A corresponding fi- 
nite number of sampling points suffices for reconstruc- 
tion. Since sampling theory has its origins and most of 
its applications in communication engineering, sampling 
theory on Riemannian manifolds has been little studied 
so far. An exception is the SU(2) group manifold, for 
which the spectral cutoff yields the much-discussed fuzzy 
sphere, see, e.g., [l^. Very interesting results were ob- 
tained by Pesenson, see, e.g., 0, who considered, in 
particular, the case of homogeneous manifolds. In jlfij . 
reconstruction works differently, however, namely by ap- 
proaching the solution iteratively in a Sobolev space set- 
ting. Useful methods should also be available from the 



field of spectral geometry, see, e.g., |17|. which studies 
the close relationship between the properties of a mani- 
fold and the spectrum of its Laplacian and, in particular, 
from the field of noncommutative geometry and the tech- 
niques based on the spectral triple, see [lSj . 

Let us now turn to the entirely new case of sampling 
in pseudo-Riemannian manifolds. We define a (N+l) di- 
mensional covariant "bandlimit" as a cutoff of the spec- 
trum of a scalar self-adjoint differential operator such 
as the Dirac or the d'Alembert operator □ (as opposed 
to, e.g., regularization through a generalized £- function). 
While these operators are self-adjoint, they are not ellip- 
tic and their spectrum needs to be cut off from above and 
below. This will lead us to a new sampling-theoretic prin- 
ciple that accounts for the Lorentz contraction of lattice 
spacings: Each temporal frequency component 4>(uj,x) of 
a field possesses its own finite spatial bandwidth and can 
be reconstructed from discrete spatial samples of corre- 
sponding density. Equivalently, each spatial mode 4>(t, k) 
possesses a finite temporal bandwidth and can be recon- 
structed from samples taken at correspondingly dense 
discrete sampling times. 

To see this, consider first the case of fiat (N+l)- 
dimensional Minkowski space-time. Fourier theory is 
applicable and cutting off the spectrum of, e.g., the 
d'Alembertian amounts to requiring \pq — p 2 \ < lp 2 , as 
is implementable in the action, e.g., through a power se- 
ries in Pq — p 2 with radius of convergence lp 2 . Thus, for 
each fixed po there is a finite bandwidth volume, B(pq), 
in the A-dimensional space of p values. The bandwidth 
volume is of ball shape if \po\ < lp 1 and is of spherical 
shell shape if \po\ > lp. Each frequency component po 
of a field possesses its own finite spatial bandwidth and 
can be reconstructed from discrete samples: 

0{po,x) =^2(f>{po,Xfi(po))G(po,x,Xfi(p )) (6) 

n 

For example, a kernel, G, for stable reconstruction from 
equidistant samples follows from the sampling theorem of 
Eqfflby viewing each of the bandwidth volumes as con- 
tained in a rectangular bandwidth box. By Landau's the- 
orem, |l3| , the precise minimum sample density for stable 
reconstruction of a temporal frequency mode, po, is given 
by the bandwidth volume B(p ) in momentum space. For 
sub-Planckian frequency modes, \po\ < lp, the maxi- 
mal bandwidth volume B(po) is the volume of the D- 
dimensional sphere with radius \Z2lp 1 . Thus, all these 
modes can be reconstructed stably from samples spaced 
at around the Planckian density. For trans-Planckian fre- 
quency modes, |po| lp , the scaling of the bandwidth 
volume depends on D, namely B(p ) = 0(|po| D ~ 2 ) for 
\po\ —>■ oo. Thus, by Landau's theorem, the minimum 
sample density for stable reconstruction as |po| —> 00 is 
decreasing for D = 1, while constant for D = 2 and in- 
creasing for D > 3. Therefore, in one and two spatial 
dimensions, Planck scale sample density suffices to sta- 
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bly reconstruct all temporal frequency modes, i.e. the 
entire field. In the curved space setting this might apply, 
e.g., to strings/string bits and to the holographic princi- 
ple for horizons respectively, see, e.g., [ljj. For D > 3 
no single sample density suffices to stably reconstruct si- 
multaneously all temporal modes. Similar to the discus- 
sion above, it is clear that each spatial mode p possesses 
a finite temporal bandwidth and that, independently of 
D, all spatial modes can be stably reconstructed from 
times series of Planck scale spacing. Generally, if the 
reconstruction is not required to be stable, significantly 
sparser sample densities may suffice. 

The analysis of Minkowski space generalizes straight- 
forwardly to static space-times. These possess coordi- 
nates in which the d'Alembertian consists of a simple 
temporal part and an elliptic self-adjoint spatial part. In 
this case, each temporal frequency mode possesses the 
equivalent of a spatial bandwidth, and vice versa. It 
should be most interesting to carry through a correspond- 
ing analysis of the sampling theory in generic space-times 
with singularities and with horizons, where space-like and 
time-like coordinates can switch their roles. Recall that 
while the reconstruction formulas of fields depend on the 
choice of coordinates, the physics of the bandwidth cut- 
off is, of course, invariant. For those pseudo-Riemannian 
manifolds which allow Wick rotation the situation re- 
duces to the simpler case of Riemannian manifolds and 
Wick rotation of the reconstruction formulas. 

Of interest to phenomenology is that, in typical infla- 
tionary scenarios, modes which are today of cosmological 
size grew from the Planck scale merely about five orders 
of magnitude before their dynamics froze upon cross- 
ing the Hubble horizon. Therefore, Planck scale physics 
could have a small but potentially measurable effect on 
the inflationary predictions for the CMB amplitude- and 
polarization spectra, see A small effect is predicted, 
in particular, if the above-mentioned generalized uncer- 
tainty relations hold, see , which would give fields the 
sampling property. So far, these models introduce Planck 
scale physics in ways that are tied to the preferred folia- 
tion of space-time into the essentially flat space-like hy- 
persurfaces defined by the cosmic microwave background 
rest frames, and these models are, therefore, breaking 
general covariance. The new approach here can be used 
to predict possible effects of a generally covariant sam- 
pling theoretic cutoff on CMB predictions. While free 
field theory suffices for this purpose, note that field the- 
oretic interaction terms, i.e., higher than second powers 
of fields (second powers occur as scalar products in the 
Hilbert space of fields) would have to be nontrivial in 
order to yield a result within the cut-off Hilbert space. 

Interestingly, regularization of field theories in curved 
space is known to induce in the action a series in the 
curvature tensor, see [2(3] . We now see that if this series 
possesses a finite radius of convergence it induces a band- 
width cutoff and sampling theorem for the metric itself, 



i.e., a cutoff on the curvature of space-time. 
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